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O , Abstract. In [Proc. Roy. Soc. London Ser. A 453 (1997), no. 1962, 1411-1443] A.S. Fokas 

04 ■ introduced a novel method for solving a large class of boundary value problems associated 

, with evolution equations. This approach relies on the construction of a so-called global 

■ relation: an integral expression that couples initial and boundary data. The global relation 
can be found by constructing a differential form dependent on some spectral parameter, 

' that is closed on the condition that a given partial differential equation is satisfied. Such 

\ a differential form is said to be fundamental [Quart. J. Mech. Appl. Math. 55 (2002), 457- 

479] . We give an algorithmic approach in constructing a fundamental fc-form associated with 
a given boundary value problem, and address issues of uniqueness. Also, we extend a result 
of Fokas and Zyskin to give an integral representation to the solution of a class of boundary 
value problems, in an arbitrary number of dimensions. We present an extended example 

■ using these results in which we construct a global relation for the linearised Navier-Stokes 
\ equations. 
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1 Introduction 

The traditional analysis associated with boundary value problems (BVPs) are the methods 
of classical transforms, of images and the use of fundamental solutions. These have been an 
incredibly successful set of tools and cover a large class of problems. They are not however, 
infallible - an example would be the heat equation on the half line with moving boundary [5] . 
A new, unified approach to studying BVPs was introduced by A.S. Fokas in [3] which allows the 
^ ' study of problems for which previous methods would fail. This method makes use of the classic 

■ integral theorems and elementary complex analysis - using which it provides: 

• existence results; 

• well posedness conditions; 

• exact solutions in terms of Fourier type integrals. 
The last of these points deserves elaboration. In the case of evolution equations on the half line 

dtq{x, t) + u {-id:,) q{x, t) = 0, {x, t) G = M+ x (0, T), 
the new method provides an exact solution which takes the integral form 

JdD+ 

where qoik), g{k) are known functions related to the initial and boundary data respectively. 
The contour dD~^ is determined by the polynomial i.v{k), and lies in Qk > 0. This solution 
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is open to asymptotic analysis by elementary means, since all the {x, t) dependence is in the 
exponential term. Also, the contours can be adjusted by use of Cauchy's theorem to give uniform 
convergenc^. Numerical implementation of these solutions is then much simpler. For example, 
the classical solution to the heat equation on the half line is given by a Fourier series, whereas 
this method provides a representation in terms of uniformly convergent Fourier type integrals 
which are easily dealt with using standard numerical techniques. 

This method of solving boundary value problems is reliant on the analysis of the global 
relation. This is an integral relation, dependent on some spectral parameter, that couples all 
the data on the boundary. To construct this, one must first recast the governing PDE into 
a divergence form, which is then open to applications of classical integral theorems. 

We treat this construction in the setting of the exterior calculus, so by means of Stokes' 
theorem we may proceed to higher dimensions arbitrarily. Within this setting, we proceed 
formally: given a boundary value problem on with associated PDE Cq = 0, we construct 
a differential form rj = ri{q,a) G A^{Q) such that 

dr/ = ^ Cq = 0, 

where a G C is a spectral parameter. Such a differential form is said to be fundamental [5]. We 
see that under the condition Cq = in Q, Stokes' theorem gives 



which is the global relation. So for each global relation, there is a corresponding fundamental 
fc-form and we use this as a basis for constructing global relations for systems of linear PDEs. We 
give a specific example by constructing the global relation for the unsteady Stokes' equations. 

We then use the construction of the fundamental /c-form to extend a result due to Fokas 
and Zyskin concerning integral representations of solutions to BVPs in two dimensions. This is 
achieved by proving their result in the context of the work in this paper, and then extending it 
to arbitrary dimension. 

We confine our attention to linear partial differential operators (PDOs) with constant coeffi- 
cients, but the majority of the analysis can be carried over into non-constant coefficients rou- 
tinely. More generally, one can work with a ring of pseudo-differential operators (see [l] for 
example), but we shall not pursue this line of work here. 

1.1 Organisation of the paper 

In Section[2]we give a formal outline of the construction of a fundamental A;-form for a given BVP. 
This involves a series of elementary lemmas, the results of which can be applied algorithmically. 
We give an example of the construction for some simple differential operators, and outline the 
use of the global relation in solving the classical problem associated with travelling waves on 
a string. A result from [5] is then extended to give an integral representation of the solution 
(assuming existence) to a class of BVPs on C M" in terms of the fundamental /c-form. We 
carry over the analysis to deal with systems of PDEs in arbitrary dimension, and conclude 
in Subsection 12.41 with an extended example in which we construct the global relation for the 
unsteady Stokes equations in 3+1 dimensions. We aim to keep notation standard throughout, 
but the reader is referred to the appendix in case of confusion. There, we also include some of 
the longer calculations that are referred to in the main body of the paper. 

^Certain conditions must be placed on the function class to which qo{k) belongs, but not overly restrictive. 
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2 The fundamental fc-form 

In this section we introduce the concept of the fundamental fc-form: this is a differential form 
associated with a given linear PDE such that the form is closed if the PDE is satisfied. Use of 
a solution to the adjoint problem then allows us to enforce that the differential form is closed 
if, and only if the PDE is solved. This then gives an equivalent problem, involving integrals 
along the boundary of the domain of our PDE ~ thus naturally incorporating the boundary data 
for the problem. We study these integral equations, or so-called global relations, and invoke 
spectral analysis to completely solve BVPs associated with systems of linear PDEs [3]. 

2.1 Formulation of the problem 

We define 7] € A"~^ (Q) for the linear PDO >C = ^ Cq9" (multi-index notation assumed) via 

dr] = {qCq - qC^)dx^ A • • • A dx", (1) 

where il. C is a simply connected domain with piecewise smooth boundary and the 
are constants. Later, we will introduce {q, q} to be such that Cq = C^q = and which gives 
a sufficient condition for ij to be closed: that is dr/ = 0. An application of Stokes' theorem then 
gives 

J an 

It will be this equation that links the given data on the boundary of our domain, to the solution 
of our problem by introducing a system of spectral parameters {aj} such that q = exp [iajx^) 
satisfies C^q = 0. These will allow us to invoke spectral analysis to solve the problem. It 
is possible to set this work up in full generality, using the theory of jet bundles to describe 
differential operators on manifolds. However, this introduces technicalities that do not aid the 
presentation of the theory that follows - so we shall not pursue it here. Essentially, the important 
property of the domain we have decided to work on is that it is cohomologically trivial: i.e. all 
closed forms are exact. 

Since r] is closed, there 30 £ A"'~^(Q) such that rj = dO, which is simply a consequence of the 
Poincare lemma and the fact Q C M" is contractible. However, it seems easier to construct 

7? = ^j{-iy+^aj{x)dx^ A • • • A dxJ A • • • A dx"", 

where the hat indicates that the particular 1-form be removed. Indeed, this formulation leaves 
us to find the {aj{x)} such that 

Ej djaj{x) = qCq - qC^, (2) 

where dj = d/dx^ . The aj{x) necessarily exist, assuming the existence of the adjoint, . This 
follows from the definition of the adjoint, since 

for any test functions (j),ip £ C^{Q,), and (•, •) is the inner product associated with L'^{Q). Our 
task now is to construct the closed form r], given C. Then, using a solution to the adjoint 
problem, which incorporates spectral parameters, the condition dr/ = is then equivalent to 
Cq = 0. 
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Definition 1. For a,/3 € Z", the anti-symmetric object [a,/3] is defined via 

[a,/3] = d'^qd''q-d(^qd"q 
and the symmetric object {a, (3} is defined by 

{a, (3} = d^qd'^q + d^qd'^q, 

where the functions q, q are suitably smooth. 

This notation will be used extensively in what follows: by proving a small collection of 
lemmas, we have means of decomposing the RHS of ([2]) into a divergence form, hence finding 
the unknown functions aj{x). The convenience of this notation is realised after observing the 
following result. 

Lemma 1. Each PDO of the form >C = ^ Cq(9" decomposes as 
C = Cq ® Cq, 

where Cq = C\ and Cq = —Cl, i.e. they constitute the self-adjoint and skew-adjoint parts of the 
operator C. 

This lemma will allow us to deal with the symmetric and anti-symmetric terms that will 
appear in the sum separately, and simplify the analysis somewhat. The proof of the lemma is 
trivial: simply split the sum into parts in which |a| is even/odd. We are now left with the task 
of decomposing the sum into a divergence form, and using our new notation we have 

{qjCq-qj[:^q)= c,{a,0}+ J] c„[a,0]. (3) 

\a\ odd lal even 

We now address the decomposition of these two terms separately. 
Lemma 2. For 1 < k < n we have the identity 

[a,P] = dk{[a - ek,f3] - [a- 2efc,/3 + e^] H 

+ (-l)(^'=+^)[a - 7fcefc, /? + (7fc - l)efe]) + {-l)"'^[a - 7^6^,/? + 7^6^] 

and the analogous result for {•, •} is 

{a, (3} = dk{{a - ek,l3} - {a - 2ek,P + e^} H 

+ (-l)(^^+^){a - 7fcefc, /3 + (7^ - l)efc}) + " TfcCfc, /? + 7^6^} 

for any 7^ < ak, where a = (ai, . . . , 0;^, . . . , a^) and the {ej} are the standard basis vectors 
on W. 

Proof. We proceed by an elementary observation 

d^qd'^q = L»fi . • • . • • D^"qD'^' • ' ' ^fc' " ' ' 

= DI{d{^ • • • I?^ • • • DtqD^- . . . D^^-' • • • D^-q) 
- L'f • • • • • • D^^qD^^ ■ ■ ■ Z)^*-^ • • • D^"q, 

where = (^d/dx^)"^. Clearly then, we have 
[a,P] = dk[a-ek,/3] - [a - e^, /3 + e^]. 
The result follows, inductively. The proof for {•, •} is similar. ■ 
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Lemma 3. For a € and 1 < k < n we have 

{a,a + ek} = dk{a,a}, 
where the are standard basis vectors on M". 
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Lemma 4 (Exchange lemma). For a S and e^ the standard basis vectors in M" we have 

Qa+ei^ hem q Qa+em+i-\ \-e2mq _ Qa+ei^ hSfeH hem+em+j ^ Qa+Sm+iA heJ^T+TH he2m+efe- 

_j_ ^Qa+eiA he^feH hem ^ Qa+Cm+iA hfim+jH he2mg:j 

— dm+j l^f}""'"^^"' hem q ga+em+i-i hS^I+TH he2m+efe g;j 

where I < j,k < m. 

Remark 1. In particular, this exchanges the partial derivatives between q and q, modulo some 
total differential, so repeated application gives 

[a + ei H h Cm, Oi + Cm+i H h e2m] = i9(terms involving q, q, dq, dq) 

since the final term, after exchanging each of the derivatives on the first term in \a+- ■ ■ , a+- ■ ■], 
cancels with the second term in [« + ••• , a + • • • ] . 

The proof to these all follow from the definitions, so for the sake of brevity they are omitted. 
We now use these lemmas to give an algorithmic method of reducing the sum in ([3]) into a di- 
vergence form. In what follows will denote the k-th component of the vector a € Z" , and 
|q| = ai + • • • + an- It suffices to prove the result for a generic term involving [a,0] or {a,0}, 
since the general result follows by linearity. 

We observe that if each component of a is even, then \a\ is even and we can use Lemma [2] 
on [a,0] exactly times for each a^, so that the final term will be [^a, ^a] which is zero by 
the anti-symmetry of [•,•]. Alternatively, we may have \a\ even if all but some even subset of 
the afc, 1 < /c < n are odd. Let this subset be of size 2m, 1 < m < ^n. We can make repeated 
use of lemma [2] to decompose the object into something of the form 

[7 + ei H h em, 7 + ^m+i H h e2m] + 9(other terms), 

2m 

where we have labelled a = 27 -|- ^ without loss of generality. This final term can be 

k=l 

decomposed using the exchange Lemma [H This concludes the proof of the existence of the 
decomposition in the case for \a\ is even. There is only one case to consider for \a\ odd: that is, 
when there is an odd sized subset of the that are themselves odd. Once again, using Lemma 
[2] repeatedly, we may decompose our object into 

{7 + ei H h em, 7 + Cm+i H h e2m+i} + <9(other terms) 

and then apply the exchange lemma on the first term, m times. Then, modulo total derivatives, 
we are left with the remainder from the exchange lemma, and the second term from {a + - ■ ■ ,a + 
■ • •}■ These two terms can be written as {/?, (3 + e^} for a suitable /? having chosen e^, and by 
Lemma [3] we are done. This construction asserts the following theorem: 

Theorem 1. Let C be a linear differential operator with constant coefficients. Given data on 
the (piece-wise smooth) boundary of a domain Q C K", 3r] £ A'^~^{Q) such that rj is closed iff 
qCq — qC'^q = 0. Moreover, rj is local and bilinear in q, q and derivatives thereof. 
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This essentially follows from the construction we outlined previously. Up to now, we have 
dealt in full generality - in practice, performing the decomposition is a simple task of using 
Lemmas EJ [3l H] repeatedly, there is very little computation involved. We proceed by way of an 
example to highlight the methods described thus far. 

Example 1. We consider the wave equation in 1+1 dimensions with C = df — d'^, on a finite 
interval - say {x, t) G 0, = (0, 1) x (0, T) with initial and boundary conditions 

q{0, t) = q{l, t) = 0, q{x, 0) = n(x), qt{x, 0) = vix). 

Then employing the methods outlined previously, we find 

= qCq -qC^q = qqu - qqtt - {qQxx - qQxx) 
= dt {qqt - qqt) - dx [qqx - QQx) , 

from which we can introduce the closed 1-form r] € A^(n), defined by 
r] = im - qqt) dx + {qq^ - qqx) dt. 

Using two solutions for^ = exp{iA;(x it t)}, we find two global relations by integrating about 
do, for the two differ eniG q 

e-'^' {qt{k, t) + ikq{k, t)) = g{k) + ikf{k) + e'''^hi{k, t) - h2{k, t), (4) 
e'''^ {qt{k, t) - ikq{k, t)) = g{k) - ikf{k) + e'^^hi{-k, t) - h2{-k, t), (5) 

where we have invoked the boundary conditions at x = and x = I. These global relations are 
valid for all /c G C, so in (l5|) we are free to take k i— > —k, and upon subtracting ^ we find 

e-'^^ {qf{k,t) +ikf{k,t)) = g'{k) + ikf{k) + sm{kl) hi{k,t). 

If we evaluate this expression for kl = nir, n d'L then we eliminate the unknown data and find 

dt{(^^'q\Kt)] = e^'^'{g\k) + ikfik)}, k = ^. 

Integrating up, we find the classical Fourier series solution to the problem. 

The previous example shows how this method works as expected for basic boundary value 
problems, but does not make much of the use of the main result so far; the next example will 
illustrate its use. 

Example 2. By means of a more involved example, we consider the differential operator C 
defined by 

^ — dxxdyydzz ~\~ Qxxdyy ~l~ Qzz- 

Equivalently, introducing the vectors Vj G via 

z;i = (2,2,2), t;2 = (2,2,0), ^73 = (0,0,2), 
we may employ the same notation introduced earlier, so that 

qCq-qC'^q=[vi,^] + [v2A + [nA- 
^See the Appendix for an outline of this computation and definitions of the notation used. 



The Fundamental A;-Form and Global Relations 



7 



We now explicitly compute the first of these terms, with repeated use of the result in Lemma [2l 
We have 

[1.1,0] = a,. [(1,2,2), 0] - [(1,2, 2), (1,0,0)] 

= dx [(1, 2, 2), 0] - dy [(1, 1, 2), (1, 0, 0)] + [(1, 1, 2), (1, 1,0)] 
= dx [(1, 2, 2), 0] - dy [(1, 1, 2), (1, 0, 0)] + [(1, 1, 1), (1, 1, 0)] , 

where the final term [(1, 1, 1), (1, 1, 1)] has been discarded, using the anti-symmetry of [•,•]. 
Making similar calculations for the last two terms, we find 

— Ox {QQxyyzz QQxyyzz ~l~ QQxyy QQxyy} dy {QxQxyzz QxQxyzz ~l~ QxQxy QxQxy} 
+ dz {QxyQxyz " QxyQxyz + QQz - QQz} ■ 

Now introducing a solution q = exp {iGjX^ ) to = 0, we find the constraint 

2 

This constraint has a rational parameterisation: introduce A G C, to find 

aia2 = ±2(A-i)"\ ctq = ±2 (A + i) . 

This can then be used in our expression for ry, from which we can again invoke Stokes' theorem 
to get an integral equation involving the spectral parameter A E C, with some additional E C 
and using similar methods to the previous example (although this case will be considerably more 
laborious), we can solve the problem given sufficient data on the boundary of our domain. 

Remark 2. Note that the construction of r] in this proof does not provide a unique divergence 
form for the differential operator - the construction is dependent on the order in which we apply 
the result of Lemma [2] to each component of a, as well as how we apply Lemma HI 

Clearly, the lack of uniqueness means there is a certain amount of ambiguity in how one should 
construct r] for a given PDE. In general, it should be chosen so that the resulting integral equation 
contains terms which are specified as much as possible by the given data on the boundary. We 
illustrate this fact by counting the number of different ways to decompose a single term in ([3]). 

m 

Without loss of generality we set a = 27 + ^ e.j for some < m < n. 

1=1 

Definition 2. Given a E , we define a path, 7'(a), to be an ordered collection of basis vectors 
such that [a, 0] (alternatively {a, 0}) is decomposed into Y) + (resp. {7 + X]j^i'7}) t>y 

applying Lemma [2] to the terms in 27 in the order V{a). 

Example 3. Let a = (2, 2, 4) = 2(1, 1, 2), so that [a, 0] can be decomposed by applying Lemma[2] 
to the terms {61,62,63,63} in some order. For instance, applying it along the path Vi{a) = 
{61,62,63,63} yields 

[a, 0] = di[ei + 262 + 4e3,0] - 52 [61 + 62 + 4e3, 61] + 93[6i + 62 + 863, 61 + 62] 

- (93[6i + 62 + 263,61 + 62 + 63]. 

Alternatively, if we apply the lemma along the path 7^2 («) = {63, 63, 61, 62} we find 
[a,0] = a3[26i + 262 + 363,0] - d^[2ei + 2e2 + 263,63] + ai[6i + 262 + 263,263] 

- d2[ei + 62 + 263,263 + 61]. 
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This example illustrates how the number of different paths for a term in the PDO will effect 
the number of different fundamental A;-forms one may construct by means of our algorithm. 
Indeed, we see that Lemma [2] can be applied in a variety of ways to a single term, the total 
number of which would be equal to the total number of distinct permutations of a generic path. 
For instance, in the previous example we could have constructed a total of 4!/2! = 12 distinct 
paths, which is exactly the number of proper permutations of the set {ei, 62, es, 63}. 

Definition 3. Given a € and a path 'P(a) we denote the number of distinct permutations 
of the elements of the path V{a) by a{a). 

Example 4. Given a = (2, 2, 5, 6) = 2(1, 1, 2, 3) + (0, 0, 1, 0), a path must consist of the elements 
{ei, 62, 63, 63, 64, 64, 64}. It follows that a{a) = 7!/(2! x 3!) = 420. 

Remark 3. Note that the length of a path will be equal to ^ (|a| — ^ of odd components in a). 

Of course, there is more to be done if some of the ak, 1 < k < n happen to be odd: we need 
to apply the exchange lemma a definite number of times. We know that we can apply Lemma [2] 
successively to reduce a generic term into either 

[7 + eiH he2m,7] (6) 

in the case of an even number of odd components of a, or into 

{7 + eiH he2m+i,7} (7) 

in the case of an odd number of odd components of a. Concentrating on ([6]) first, we see there 
are (^^) ways to manipulate this object, by means of Lemma [21 into the form 

[7 + efci + ■ ■ ■ + efc^,7 + ekm+i + ^ ^k^J, 

where the {e^j} are some permutation of the {ei}. Now in this form, we have mx m choices for 
the first application of the exchange lemma, then (m — 1) x (m — 1) for the next, etc. So the 
total number of ways of applying the exchange lemma to this term is (m!)^. And so the total 
number of ways of decomposing ^ into divergence form is 

(2m)! 

X ml X ml = (2m)!. 



ml X ml 

Similarly, there are {^"^^) ways to manipulate ([7]) into the form 

{7 + efci + • • • + efc„, 7 + efc,„^i + • • • + e^^^^ J 

by means of Lemma [2l We then have (m + 1) x m choices for the first application of the 
exchange lemma, followed by m x (m — 1) choices etc. It follows that the total number of ways 
of decomposing ([7|) into a total divergence is 

(2m + 1)' 

--^ . X (m + 1)! X m! = (2m + 1)1. 

(m + 1)! X m! 

In both cases then, the total number of different decompositions of the terms ([6]) and ([7]) is 
simply of odd terms in a)l. We can now state the total number of fundamental /c-forms that 
can be constructed using out method. 

Theorem 2. Given a linear PDO with constant coefficients C = X]qCq,(?", our method can 
construct N(C) different fundamental k-forms, where 

N{C) = l[OJa{a) 

a 

and Oa denotes the number of odd components of a G 
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The proof follows from the arguments previous based on a single term. The total number is 
found by simply taking the product over all the terms in the PDO. From this result it immediately 
follows that in Example [H the fc-form we constructed was one of 

N{C) = 1x1 = 1 

for C = dtt — dxx- So in this case, it was in fact the unique fundamental form for the BVP. 
However, for the higher order operator in Example [2l where £ = dxxdyydzz + dxxdyy + dzz-, the 
total number of distinct A;-forms we could construct is 

N{C) = 3! X 2! X 1! = 12, 

any of which could be used to solve the BVP, given sufficient data on dVt. 

Example 5. We consider bi-harmonic functions on W', i.e. classical solutions to A^g = 0. In 
Cartesian coordinates then, the relevant differential operator is given by 

^ dxXXX ~l~ 9yyyy + QzzZZ ~\~ ^Oxxyy ~l~ ^Oyyzz ~\~ ^OzzXX- 

Then employing the notation introduced earlier, we seek to decompose 

[^;i,0] + [V2,0] + [t;3,0] + 2[^;4,0] +2[7;5,0] +2[7;6,0], 
where the {vj} are defined by 





= (4,0,0), 


Vi 


= (2,2,0) 


V2 


= (0,4,0), 


Vb 


= (0,2,2) 


V3 


= (0,0,4), 


V6 


= (2,0,2) 



From which, we can read off the number of distinct fundamental /c-forms produced by our 
method to be 

N{C) = 1! X 1! X 1! X 2! X 2! X 2! = 8. 

Below is an example of one such decomposition that can be used to construct the fundamental 
form r], where the {aj} are defined so that J2j (^^ adjoint solution q = exp{iajX^) 

solves A^g = 0) 

= dx {e'"'''' [qxxx + icrfq - idiqxx - crfqx + 2qxyy + iaialq + 2qxzz + 2iaialq]^ 

+ dy i^e^^j"^^ {qyyy + ialq - ia2qyy - ajqy - 2iaiqxy - 2aia2qx + 2qyzz + 2iaialq}^ 
+ dz {e""'"^' {qzzz + io-lq - ia^qzz - crjqz - 2iaiqxz - 2aia2,qx - 2ia2qyz - 2fT2fT3(?y}^ , 

where {x,y,z) = {x^ , x'^ , x^^) . In constructing this equation, we simply made repeated use of 
Lemma [2]- there is very little computation involved. 

We must now address the question of whether N{C) is well defined - that is to say: is N{C) 
a coordinate independent object? Suppose Si is the set of fundamental /c-forms associated 
with C on the domain Oi, and let p : fii ^ ^2 be a diffeomorphism, with inverse (j). Then the 
question becomes: is N{C) on Qi as on il.2- This is partially addressed in the following lemma. 

Lemma 5. Each fundamental k-form in ill induces a corresponding fundamental k-form in 0.2 
and vice versa. 
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Proof. Take € S2, i.e. rj is fundamental on i}2- Then the pull-back p* induces a fundamental 
form on Oi. Indeed, setting 6 = p*rj, we find 



since the pull-back commutes with the exterior derivative. We conclude that for each element 
of S2 there is a corresponding one in Si. Using the same argument with the pull-back (p*, 



The reason this doesn't fully answer our question, is because of the following. Let t], 9 be 
two fundamental A;- forms on Q constructed via our method and consider (f>*r] and (p*9. These 
are fundamental on but are not necessarily of the form constructed via our method. That 

is to say, applying our method to C after applying the diffeomorphism to fi, wont necessarily 
yield the same set of fundamental A;-forms as it would if we computed the induced fundamental 
/c-forms via the appropriate pull-back. 

Now although we have shown that many fundamental forms may be constructed, they are not 
necessarily independent. We use the fact that there is no cohomology on the class of domains 
we are working on, so any closed form is exact. Now consider two fundamental /c-forms rj, 9 that 
differ by some closed (hence exact) form dip. These two forms yield the same global relation, 
indeed 



since d^ = 0. We should confine our attention then, to the equivalence class of fundamental 
forms in which r^i ~ r/2 rji — r]2 = d9, for some 9 € A^~^{Q). Labelling the set of exact 
differential forms on Q by 11, and the set of fundamental fc-forms on 17 by S, we have the 
following result. 

Lemma 6. There is a unique fundamental k-form in S/II. 

Proof. Assume that S/IT contains m > 1 fundamental /c-forms, and we consider the difference 
of two such A;-forms: ip = rji — rjj, i ^ j. Note that dtp = regardless of whether qCq = qC^q. 
Since we are working on a cohomologically trivial space, it follows that 3cj G A^'~^(i7) such 
that ■0 = dw. However, this implies that rji and r]j are equivalent in S/II, contradicting our 
hypothesis. This is true for each pair so we may conclude that S/II contains one unique 

element. ■ 

The consequences of this fact are as follows. Even though we have shown that we can 
construct several different fundamental A;-forms for a given PDE, they are all equivalent up 
to some closed form, the closeness of which does not depend on the condition qCq = qC^q. 
Therefore, any of the fundamental /c-forms we construct will generate the same global relation. 
In practice then, we should choose a particular decomposition that yields a fundamental /c-form 
whose coefficients contain derivatives of q that match as closely as possible, the prescribed data 
on the boundary. This will then lead to a global relation that can be dealt with efficiently, by 
using the appropriate data on the boundary. 

2.2 Integral representations 

Here we use the fundamental fe-form to give an explicit solution to a EVP associated with the 
polynomial differential operator C = £{81,82, ■■■ ,dn)- We expand on the observation given 
below: 



d^ = d {p*r]) = p*dr] = 



completes the proof. 
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Lemma 7 (Fokas & Zyskin). Given a polynomial differential operator p{di, 82), and rj € A-'^(i7) 
such that 

dri = p(i/ci, iA;2)g(xi, X2)e"*''i''^"'''2^2 ^^^^ 
where p{di,d2)q = in il., then if q exists it can be represented in the form 

A proof of this is found in [5], it essentially follows from the distributional relation 

We can now extend this result to arbitrary dimension as follows. We construct an fundamental 
k-form r] using the lemmas of the previous section, but we do not solve the adjoint problem. 
Instead, we introduce q = exp{—ikiX^) from which it follows that 

C^q = Ci-di, -d2, . . . , -5„)e-^^'"' = C{ik,,ik2, . . .,ikn)e-'^'''\ 
Now using equation ([T]), we find that 

dr/ = -C{iki,ik2, ■ ■ ■ , ikn)q{xi, ... , x„)e~*^'''' dx^ A • • • A dx"' 
assuming Cq = in Q. This then gives the integral representation for q as 

' ^ dk,...dkj y^dn. 



(27r)" yjjn Jq^ C{iki,...,ikn) 

We summarise this result in the following. 

Theorem 3. Given a boundary value problem in Q, with associated polynomial differential 
operator C = C{di, . . . ,dn), on the assumption that a solution exists, it can be represented by 

-1 /■ f e*'=»^'7?(y,/c) 

= /o \r, / dfcl • • • dkr, 



(2vr)" jQQ£{iki,...,ikn)' 
where y £ dQ and rj is a fundamental k-form for C associated with the adjoint solution e~**'»^\ 

Remark 4. Note that this representation is in accordance with the Ehrenpreis principle from 
the theory of complex analysis of several variables. See for instance [2]. 

The term 'adjoint solution' is a misnomer, q does not solve the adjoint problem. It is simply 
the function we choose for q in the decomposition of ([3]). Note that the non-uniqueness of rj 
discussed earlier, does not affect this result. Indeed, the addition of some exact form to r]{x, k) 
will not contribute since Jg^ d'0(x, k) = 0. 

2.3 Systems of PDEs 

Here we extend the results proven in the last section, to deal with systems of coupled PDEs. The 
vector (p = {4'i,(p2j ■ ■ ■ :4'n) will be denoted by 4>i and we shall assume summation convention 
throughout. The differential operator £ will be expressed component wise, so that we denote 
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£ = Cij. Studying linear systems of the form C[ 
linear set of PDEs 



then translated to solving the coupled, 



0. 



J~'n\ ' ' ' ^nn 

From this, we construct the formal adjoint, defined by 
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nl 



.'-Ir 



We then decompose an analogous object to that in the previous section, given by 



r 



Using the definition of the adjoint, and changing the order of summation in the second term, 
we see this is equivalent to decomposing 



Now we may apply the results of the previous section to each of the Tn? terms in this sum. 
Since we will want the closeness of the fundamental fc-form to be equivalent to £[i?!)] = 0, we 
must choose the solution to the adjoint problem so that each of the scalar fields 0j are linearly 
independent functions. With this condition, we see that we must have C 



0, i = 1, 



,n. 



for the fundamental /c-form to be closed. This is exactly the condition £[(/>] = 0. 
2.4 Unsteady Stokes equations 

In this section we construct the fundamental 3-form associated with the unsteady, unforced 
Stokes equations, which are a direct linearisation of the classical Navier-Stokes equations. We 
use a separable to solution to the adjoint problem suitable for geometries which lend themselves 
to Cartensians. We deal with a divergence free velocity field u(x, t) = (ui, U2, U3) with and some 
scalar field p(x, t), such that 



dtui — uAui + dxP = 0, 
dtU2 - uAu2 + dyp = 0, 
dtU3 - vAu^ + dzP = 0, 
dxUi + dyU2 + dzUs = 0, 



(8) 
(9) 
(10) 
(11) 



where v > is constant and A is the Laplacian on M^. We now introduce the 4- vector t) = 
{(j)i,(j)2,<p3,4'i) where <j)4 = p(x, t) and (j)i{x,t) = Ui{x,t) for i = 1,2,3. With this notation, 
equations ([8|)-(|lip take on the form £.(j) = 0, where 



def 



dt-uA 








dx 





dt-uA 





dy 








dt-i^A 


d. 


dx 


dy 








We now employ the methods introduced in the last section to analyse this coupled system of 
linear PDEs, and see how the fundamental 3-form associated with the problem can give insights 
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into the associated IBVP. Note first 











-dx 





-dt - z/A 





-dy 








-dt - uA 


-d. 


-dx 


-dy 


-d. 






We first solve the associated adjoint problem, defined by >C^<^ = 0, incorporating a spectral 
parameter. In doing this, we are able to integrate this closed form over our domain from which 
Stokes' theorem allows us to realise a constraining equation involving Fourier-type integrals of 
data on dQ. 

Lemma 8. A solution to the adjoint problem of the required form is given by 

where k = k(,^) is an isotropic vector in C^, = (■^1,^2) is a 2-spinor and G C. 

This result may be checked routinely. The appearance of the spinor has an important conse- 
quence, as we shall now see. 

Lemma 9. There are at least 2 global relations associated with the IBVP for the unsteady Stokes 
equations. 

Proof. This is essentially a consequence of the fact Spin(3) = SU (2) is (defined as) a double 
cover of SO (3). Rotating k by 27r about some axis must leave it invariant so still provide 
a solution to the adjoint problem, but the associated spinor transforms as ^ 1-^ for U € SU (2), 
with [/ / I. ■ 

We now decompose (i8l)- (|lip into a divergence form, from which we may construct the associated 
fundamental form. We omit the details of this, but it follows from Lemmas [21 [3l [H that the 
unsteady stokes equations imply 



dp 
di 



+ Div J = 0, 



where 'density' p and the components of the associated 'flux' J = (J^, J^, J^), are given by 

p = ^^ui + (i>2U2 + (t>3U3, 

= (pip + 4)ip + v{uidx4>i - <PidxUi + U2dx(t>2 - <p2dxU2 + U3dx(t>3 - (psdxUs) , 
= ^2P + 4>2P + v[uidy4)i - (pidyUi + U2dyci)2 - 4>2dyU2 + u^dy^^i - (psdyUs) , 

= (i>3P + 4>3P + T^{uidz(i)i - CpidzUi + U2d;,(f>2 - ^2dzU2 + ^iSf^z^S " hdxUs) , 

where p = p{x,t), Ui = Ui(x,t) etc. We now construct the fundamental 3-form associated with 
the problem, in which we encompass our previous solution to the adjoint problem. That is, we 
conclude that if C(j) = then the following differential form is closed 

r] = pdx Ady Adz - J^dy Adz Adt + J^dx Adz Adt- J^dx A dy A dt. 

And so is fundamental. We note that the terms in the above differential form contain the 
(^1, ^2) '^3)) which serve as spectral parameters. The global relation is then formed by integrating 
this form over using data prescribed on the boundary. This results in an equation involving 
Fourier type integrals, which lend themselves to analysis as seen in [5]. 
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Appendix A: Notation 

We will use ft to denote a simply connected subset of M" with piecewise smooth boundary dQ. 
L^(r2) will denote the Hilbert space of square integrable functions w.r.t. the Lebesgue measure 
on Q, A^(r2) denotes the vector space of differential A;-forms on $7, d : A^(r2) A^^^(r2) is the 
exterior derivative. When working with a fundamental fc-form, k = n — 1, i.e. one less than the 
dimension of the domain we are working on. We assume Stokes' theorem: on a n-manifold M 
with boundary dM, we have 




for T] G A"~^(A^) with compact support. The space of test functions on Q will be denoted C^{Q), 
meaning the space of smooth functions on Q with compact support. C will denote a linear 
differential operator, with associated formal adjoint with respect to the inner product 
on L^(r2). The Greek letters a,/3... € will be used in mutli-index notation. The letters 
k will denote numbers ranging from 1 to n. Summation convention will be used. 



Appendix B: Wave equation in 1 + 1 dimensions 

In Example [1] we construct the 1-form r] € A'^(r2), given by 
1] = imt - mt) dx + {qqx - qqx) dt. 

Letting q solve the adjoint problem, so g' = exp[iA;(x — t)], enforces the closeness of is equivalent 
to q solving = 0. We integrate this differential form on dVt^ and it follows from Stokes' 
theorem that 

f e^*^(^-*) [{qt + ikq) dx + {q^ - ikq) dt] = 0. 
JdQ, 

Now we carry out the integration, using the following notation for the initial data 

q{k,t)= r e^'=^g(x, t) dx, g{k)= e*^'^g(x, 0) dx, />) = f e'^''qtix,0) dx 
Jo Jo Jo 

and the following notation for the relevant terms on the boundary are 

hi{u;,t)= f e-'^^qx{Ut)dT, h^ioo^t) = f e-'^-q,{0,t) dr. 
Jo Jo 

We also introduce the notation f^{k) to be defined by 

f^{k) = [ sm{kx)f{x)dx. 
Jo 

With these definitions, two global relations in Example [T] follow. 
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